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Abstract 

Both the right and left eigenfunctions and eigenvalues of the linearized homogeneous Boltzmann 
equation for inelastic Maxwell molecules corresponding to the hydrodynamic modes are calculated. 
Also, some non-hydrodynamic modes are identified. It is shown that below a critical value of the 
parameter characterizing the inelasticity, one of the kinetic modes decays slower than one of the 
hydrodynamic ones. As a consequence, a closed hydrodynamic description does not exist in that 
regime. Some implications of this behavior on the formally computed Navier-Stokes transport 
coefficients are discussed. 
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I. INTRODUCTION 



Granular gases provide an appropriate context in which a number of fundamental issues 
related with non-equilibrium systems can be addressed. Primary among these is the existence 
of a hydrodynamic description, i.e. of closed, deterministic equations for the hydrodynamic 
fields, identified as the number density n, the flow velocity u, and the (granular) temperature 
T, in the case of one-component systems. A theoretical framework for deriving macroscopic 
equations from the underlying microscopic dynamics and investigating its range of validity 
is kinetic theory. In the low density limit, the Boltzmann equation for smooth and inelastic 
hard spheres or disks has been used to derive hydrodynamic equations for granular gases 



since many years ago [lH3|. Presently, explicit expressions for the transport coefficients 
appearing in the analogous of the Navier-Stokes equations are available, and they have 
been confirmed via direct Monte-Carlo simulations 0, 0]. Nevertheless, the methods used 
are formal and do not address either the existence or the context of the hydrodynamic 
description itself, although internal mathematical consistency is accomplished at the level 
of the Navier-Stokes approximation, i.e. first order in the gradients of the hydrodynamic 
fluxes. For elastic hard spheres or disks, the problem has been satisfactorily solved by 
analyzing the spectrum of the linearized Boltzmann operator jg]. On the other hand, for 
the inelastic Boltzmann equation (IBE), although the hydrodynamic eigenfunctions and 
eigenvalues have been identified in the limit of large wave vectors {?], Isj], almost nothing 
is known about the kinetic, non-hydro dynamic part of the spectrum. Consequently, the 
dominance of the hydrodynamic modes has not been established for the IBE, although it 
has been proven for some single-time relaxation model equations js]. 

Here the validity of a hydrodynamic description for granular gases will be addressed using 
a simplified IBE. Specifically, the inelastic Maxwell model (IMM) will be employed 9Hl2|. 
This kinetic equation is obtained from the IBE for hard spheres of disks by replacing the 
velocity dependent collision rate by an effective average value proportional to the thermal 
velocity. Although other choices have also been considered jsl, [lo|, here the effective collision 
rate will be assumed to be also independent of the collision angle. Moreover, attention 
will be restricted to the modes at asymptotically long wave lengths, i.e. to perturbations 
occurring also in homogeneous systems. This will suffice to establish that the hydrodynamic 
spectrum is not isolated from the rest at strong inelasticity, contrary to what is required for 
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the existence of hydrodynamics. 

Let f(vi,t) be the one-particle distribution function of particles of mass m with velocity 
v% at time t. The homogenous Boltzmann equation for the the IMM in d-dimensions reads 

a t f(v 1 ,t) = j M [v 1 ,t\f,f], (i) 

J M [ Vl ,t\f,g} ee J da J dv 2 (a-% 1 - l) f( Vl ,t)g(v 2 ,t). (2) 

Here v(t) is an effective collision frequency. Its explicit form will not be needed, being enough 
to know that v oc nT 1 ^ 2 . Moreover, Vt^ = 2-n d l 2 /T(d/2) is the total solid angle element, da 
is the solid angle element around the direction of the unit vector a, and b~ l is an operator 
changing all the velocities v± and v 2 to its right into their precollisional values given by 
b^Vip = i?i2 = v i -F (1 + ot)vi2 ■ aa/2a, with v± 2 = V\ — v 2 . The coefficient of normal 
restitution a is defined in the interval < a < 1. The hydrodynamic fields are defined in 
terms of f(vi,t) in the usual way. Then, by taking velocity moments in Eq. (CQ) it is seen 
that 

dn du dT J _. Nm 

a = a=°- 8F = -«'> r - < 3 ) 

with ( being the cooling rate given by 

C(t) = ^(l -«>(*). (4) 

As the IBE for hard spheres and disks, Eq. ([T|) has a similarity solution describing the 
homogeneous cooling state (HCS), 



f H (vi,t) = n H v d (t)x(c 1 ), v (t) = 



2T H (t) 



1/2 



(5) 



m 

where x( c i) is an isotropic function of the scaled velocity c x = Vi/v (t). Substitution of Eq. 
(ED into Eq. using Eq. Q3J yields 

2d^ ' f ClX ^ Cl ^ = J *^ c il^x]- (6) 

In the above expression, ( = Ca{t)/vjj{t) and Jm is the dimensionless version of the collision 
operator Jm- The index H is used to characterize quantities computed in the HCS. 

Next, Eq. ([T|) is particularized for small (homogeneous) perturbations of the HCS. A 
distribution 5x is defined by f(vi,t) = riHVQ d (t)[x(ci) + 5x(ci,t)], with the velocities scaled 
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relative to vo(t). Moreover, the time is expressed in terms of r given by dr = un(t)dt. 
Retaining terms up through linear order in 5\ it is obtained 

d T 5 X (c u r) = A(ci)<Jx(ci,r), (7) 
A(ci)5x(ct,r) = Jm[c 1 \x,5x} + Jm[ci\5x,x]-^-J^-[ci5x(c 1 ,t)} 

' [cix(ci)] / dC2 (t ~~ x ) 6x (C2 ' r) • (8) 

Linearization of the balance equations ([3]) around the HCS leads to 

dp du ( 89 7 ~ 

with p = (n—nH)/n>H, u> = u/vo(i), and = (T — Th)/Th- The spectrum of the above equa- 
tions is given by three points, Ai = 0, A 2 = C/2, and A 3 = — C/2, A 2 being <i-fold degenerate. 
Their perturbation for finite wavevectors (gradients) in the context of the inhomogeneous 
IMM defines the hydrodynamic modes more generally. These eigenvalues are the same as 



for the IBE for hard spheres or disks 
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In the elastic limit a — > 1, the well known result 
of a (d + 2)-fold degenerate point at zero eigenvalue is recovered. 

The operator A is expected to have the corresponding hydrodynamic modes, the remain- 
ing part of the spectrum being referred to as the kinetic modes. Of course, this terminology 
does not preclude any mathematical difference or separation between both parts of the 
spectrum. Then, the eigenproblem 

A( Cl )&( Cl ) = \Ma) (io) 

is considered. First, attention is focussed on the hydrodynamic part of the spectrum that can 
be determined as follows. Define the function F(cx, p, u}, 7) = px{yCi), with C\ — c% — u>. 
Using Eq. (JBJ) it follows that 

^ i -(C 1 F)= 1 d JM[ci\F,F]. (11) 

Differentiating this equation with respect to p, u>, and 7, taking afterwards the limit p = 
7 = 1, u) = 0, it is straightforward to show that A has the hydrodynamic eigenvalues Aj, 
% = 1,2, 3, given above, with the eigenfunctions 

dx(ci) 



fi(ci) = (d+i) x ( Cl ) + Cl 



dci 



6(ci) = -^, &(Pi) = ~ ■ [c lX ( Cl )] . (12) 



As for the hydro dynamic equations, the eigenvalue A2 is <i-fold degenerate. The eigenfunc- 
tions £j here are the same functionals of x as f° r the IBE [7] . Solutions to Eq. (J7J) are sought 
in a Hilbert space defined by the scalar product 

(g\h) = f dcx-\c)g\c)h{c), (13) 

with the dagger denoting complex conjugate. The hydrodynamic eigenfunctions span a d + 2 
dimensional subspace of the Hilbert space. They are not orthogonal, as a manifestation of 
A being non Hermitian. Then, the left hand eigenproblem 

A+( Cl )^(ci) = A^(Q), (14) 

has to be considered. In this equation, A + is the adjoint of A, defined through (g\Ahy = 
(h\A + g), 

d 1 2 /* /* 

A + (ci)0( Cl ) = / dc 2 I do- X (ciMc 2 )(6 CT -l) [x~\<*)g{ci) + X~\<%)9fa)] 

+|X(ci) ( —j- ~ 1 1 / rfc 2 x _1 (c2)^(c 2 )— [C2X(C 2 )] 

+~x(ci)ci--^[ X - 1 (c 1 )g(c 1 )], (15) 

b a being the operator inverse of b^ 1 . Although there seems to be no simple relationship 
between the functions £j and £ i5 the hydrodynamic part of the spectrum of A + has been 
found by simple inspection. The functions 

ei(ci) = x( Cl ), Ua) = c lX (c 1 ), U*)=(j + ^x(<*), (16) 



are solutions of Eq. ( jl4j) . corresponding to the eigenvalues Aj = Aj, i — 1,2,3, respectively. 
Again, the eigenvalue £/2 is (i-fold degenerate. Moreover, it is easily verified that 

= (17) 



for the hydrodynamic eigenfunctions given by Eqs. (fl2l) and ( 1161) . 

It is also possible to get some information about the remaining, kinetic part of the spec- 
trum of A + . Here attention will be restricted to two eigenfunctions which have a clear 
physical interpretation. A direct evaluation gives 

A + (c 1 )c lx c ly x(c 1 ) = + afci x ciyx{ci). (18) 



It follows that £ 4 (ci) = C\ x c\ y x{c\) is an eigenfunction of A + with eigenvalue A4 = —(1 + 
a) 2 /4. Note that £ 4 is proportional to the dynamical variable whose velocity average provides 
the pressure tensor p]. It is important to check whether the kinetic eigenvalue A4 is actually 
separated from the hydrodynamic part of the spectrum, as required for the validity of the 
hydrodynamic description. Note that if it is assumed that any function of the Hilbert 
space can be expanded in terms of the eigenfunctions of A, then the left eigenvalues are 
also right eigenvalues, i.e. there is a solution of Eq. ([9]) with A4 = A4. Consequently, the 
general solution of Eq. ([7]) contains a contribution proportional to exp(A4r)^4(ci), where 
£4(01) is the (unknown) eigenfunction of A with eigenvalue A4 = A 4 , or a linear combination 
of eigenfunctions in the case of degeneracy. The validity of hydrodynamics requires that 
the above contribution decay faster than all the hydrodynamic modes, i.e. it must be 
I A4 1 > I A3 1 = C/2. Taking into account the exact result given by Eq. ([5]), this condition is 
equivalent to a > a.\ with cti = — (d — 2) /(3d + 2). For d = 3 this expression leads to an 
unphysical negative value of oti, while for d — 2 it is a\ — 0. 

The above result suggests to check whether there is another left eigenfunction proportional 
to the dynamical variable providing the heat flux. One reason to consider these fluxes 
as candidates to be eigenfunctions of A + is that they are known to be orthogonal to the 

n 

hydrodynamic modes £j [7|. Explicit evaluation gives that 
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Z 5 (c 1 )=lcl- — )c lxX (c 1 ) (19) 
is a solution of Eq. f[T4"j) with 

The condition for the validity of hydrodynamics following from the existence of this kinetic 
mode is that a > a 2 = (4 — d)/3d. Therefore, there is a finite value of the coefficient of 
normal restitution below which there is no time scale separation between the hydrodynamic 
and the kinetic parts of the distribution function. For d — 2 it is cxi — 1/3 and for d = 3 
a.2 = 1/9. In Fig. [I] the two kinetic eigenvalues A 4 and A 5 are compared with the slowest 
hydrodynamic mode A3 for d = 3. 

The lack of dominance of the hydrodynamic modes over the mode related with the heat 
flux for long times, reflects itself in the calculation of the Navier-Stokes heat conductivity of 
an inelastic gas of Maxwell molecules. As already mentioned, the form of the hydrodynamic 
spectrum for IMM reported here is the same as for the IBE for hard spheres. The only 
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FIG. 1: The two kinetic eigenvalues A4 (dotted lines) and A5 (dashed lines), and the slowest 
hydrodynamic mode A3 (solid lines) versus the coefficient of normal restitution a for a system of 
inelastic Maxwell molecules. Note that the three modes are negative. 

difference is in the explicit expression of the distribution function of the HCS and in the value 
of the cooling rate. As a consequence, the formal expressions for the transport coefficients 
derived in [J [3] from the IBE, can be directly applied to the case of Maxwell molecules. 
Then, the time-dependent heat conductivity k is given by 



The existence of hydrodynamics to Navier-Stokes order requires that the correlation function 




(21) 



x 




(22) 



k(t) decay to zero for r ^> 1. This clearly implies that |As| > £/2, i.e. a > a<i- Consistently 
with this analysis, it is found that the heat conductivity for inelastic Maxwell molecules 
obtained by the Chapman- Enskog procedure also diverges for a — > «2 

To put the results reported here in a proper context some comments seem appropriate, 
(i) It has been proven that the solution of the Boltzmann equation for IMM tends to the 



HCS for arbitrary initial conditions 



151 ] . The lack of time scale separation discussed here 



does not contradict this general property, since the obtained kinetic eigenvalues are negative, 
therefore leading to decaying in time contributions, (ii) It is worth to stress that the failure 
of hydrodynamics in IMM should not be understood as limited to the Navier-Stokes approx- 
imation, but to any closed description of the system in terms of the hydrodynamic fields, 
(iii) Of course, the value of a for which the scale separation actually fails is not known, but 
only a lower bound has been determined. It seems sensible to expect higher limiting values 
of a associated to eigenfunctions involving higher velocity powers, (iv) On the other hand, 
there is no reason to expect a similar behavior for the IBE for hard spheres or disks. On the 
contrary, there are some indications that this is not the case. The expressions derived for 
the Navier-Stokes transport coefficients are regular functions of a for < a < 1 [l^, and 
the two kinetic modes considered here are not left eigenfunctions of the linearized inelastic 
Boltzmann operator. Actually, other relevant deep differences between the IBE and the 
IMM have been found, as for instance the asymptotic decay of the distribution function of 
the HCS for large velocities HQ- 
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